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1.  INTRODUCTION 


The  MSX  experiment  is  scheduled  to  fly  a  telescope  at  900  km  altitude. 
This  telescope  is  to  be  eryogenically  cooled  by  liquid  1*2'  teaching  a 
temperature  in  the  vicinity  of  20  K.  The  purpose  of  this  report  is  to 
provide  estimates  of  contamination  with  particular  emphasis  on  the  problems 
involved  when  the  telescope  and  its  mirror  are  facing  into  the  ram 
direction.  Two  major  mechanisms  for  contamination  are  considered: 
deposition  of  OH  formed  by  the  reaction  of  outgassing  H2  coolant  with 
atmospheric  O,  and  the  formation  of  OH  on  the  mirror  caused  by  surface 
reaction  of  atmospheric  O  and  H. 


2.  DEPOSITION  OF  OH  FORMED  VIA  GAS  PHASE  REACTION  OF  0  AND  H2 

As  ths  liquid  H2  coolant  is  ussd  up,  it  is  necessary  to  vent  the 
gaseous  H2  that  is  formed.  It  is  estimated  that  0.001  grams  per  second  of 
H2  will  be  released.  At  the  900  km  altitude,  roost  of  the  H2  will  travel 
quite  far  (on  the  order  of  100  km)  before  experiencing  a  collisional 
interaction  with  an  atmospheric  species.  The  chance  of  any  product  of  9uch 
a  collision  making  it  back  to  the  experimental  platform  is  negligible. 
However,  there  are  some  collisions  which  occur  ouch  closer  to  the  vehicle. 
If  such  a  nearby  collision  involves  atomic  oxygen,  then  a  reaction  is 
possible  forming  OH  as  a  product,  and  OH  from  such  nearby  collisions  may 
find  its  way  back  to  the  telescope.  It  is  the  purpose  of  this  section  to 
estimate  the  magnitude  of  OH  deposition  from  this  mechanism. 


2.1  Undisturbed  Ha  Flow  Field 

As  the  H2  is  released  from  its  sonic  orifice,  it  will  experience  a  few 
collisions  among  itself,  and  then  adopt  an  essentially  radic 1  flow  field 
within  a  few  exit  diameters  of  the  orifice.  Since  this  distance  is  small 
compared  to  the  other  length  scales  to  be  discussed  subsequently,  it  is 
appropriate  to  replace  the  orifice  by  a  point  source  of  exhaust  molecules 
traveling  at  their  thermodynamic  limiting  speed  (i.e.,  with  zero 
translational  temperature)  with  a  number  density  distribution  (before  any 
collisions  with  atmospheric  species)  given  by 

ru,  *  -%f(9)  »  (1) 

"2  r2 

where  B  gives  the  axial  number  density  decay  and  e  represents  the  angle 
from  the  orifice  normal.  The  particular  form  for  f(e)  that  is  used  is  an 
asymptotic  form  of  that  proposed  by  Brook, *  namely 

f (e)  -  exp{-\2 ( 1-cos (e) ]>  ,  (2) 


where 


(3) 


(5) 


(6) 


In  the  above  relations,  u0,  MQ/  nQ  and  Ag  denote  the  exit  velocity,  Mach 
number  (unity,  in  this  case),  number  density  and  area,  respectively,  and  um 
and  Y  are  the  thermodynamic  limiting  speed  and  the  exhaust  ratio  of 
specific  heats.  Several  other  forma  for  f(e)  have  been  proposed,^'3,4'5 
but  the  Brook  model  was  selected  on  the  basis  of  excellent  comparisons  with 
a  limited  number  of  method-of-characteristics  solutions.  There  has  yet  to 
be  a  definitive  comparison  of  the  various  core  flow  models  over  a  wide 
range  of  nozzle  exit  conditions. 

It  should  be  noted  that  at  large  angles  from  the  thrust  axis 
(typically  beyond  90°)  the  plume  consists  of  the  slower,  more  easily 
turned,  flow  from  the  nozzle  boundary  layer.  The  above  model  can  be 
expected  to  under  predict  the  density  and  over  predict  the  velocity  in  this 
regime.  This  is  an  unimportant  region  for  the  problem  at  hand. 

Application  of  the  above  formalism  to  the  H2  outgassing  requires 
estimating  a  temperature  at  the  orifice,  TQ.  A  temperature  of  100  X  was 
chosen  (the  results  are  not  greatly  sensitive  to  this  choice),  yielding  the 
following  results: 


uQ  *  -VYRT0/m0  •  7.6  x  104  cm/s 


(7) 


where  R  is  the  universal  gas  constant  and  mQ  is  the  molecular  weight  of 
hydrogen.  The  release  rate  of  hydrogen  molecules,  n,  is  given  by 
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*  mN  in 

n  *  u_n_A0  «  —  *  3.0  x  10^  molecules/s  ,  (8) 

u  mo 

where  m  ia  the  mass  release  rate  (0.001  grams/s),  and  N  is  Avogadro's 
number.  Substituting  Eqs.  (7)  and  (8)  into  the  earlier  equations  gives 

Ujjj  ■  1.9  x  IQ®  cm/s  ,  (9) 


B  ■  8.6  x  10^4  molecules/cm 


and 

\2  »  3.33 


(10) 


(11) 


2 . 2  Reactive  Cross  Section  for  0  ±  H3  -*  OH  +  H 

The  highest  temperature  data  that  could  be  found  for  this  reaction  is 
that  of  Roth  and  Just6,  who  measured  a  rate  constant  of 

kr  *  (3.80±0.27)xl0”10exp(-^Y^>  '  (12> 

in  the  range  of  1900-2800  K.  Given  such  a  rate  constant  with  an  Arrhenius 
temperature  dependence,  it  is  possible  to  deduce  a  unique  cross  section  as 
a  function  of  relative  collision  velocity  which  produces  that  form.7  whan 
this  expression  is  evaluated  at  a  relative  velocity  of  9.4  km/s  (7.5  km/s 
orbital  plus  1.9  km/s  Hj ) ,  it  gives  a  reactive  cross  section,  <r*,  given  by 

t t *  ■  2.1xl0"16  cm2  .  (13) 

It  should  be  noted  that  the  energy  associated  with  a  collision  at  9.4  km/s 
is  considerably  greater  than  the  2800  K  upper  limit  for  the  rate  constant 
measurement,  and  the  deduced  c*  therefore  represents  a  considerable 
extrapolation  of  the  data.  Lacking  a  better  value,  it  must  be  used  with 
caution. 
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poet-collision  relative  velocity  vector  in  terms  of  Jtnown  pre-collision 
quantities,  viz: 


*12Vii  '  2AE 


This  relation  can  then  be  substituted  into  Eq.  (14)  to  give 


_®4 _  - 

+  roA^fSrf 


where  irf  is  a  unit  vector  giving  the  direction  of  the  post-collision 
relative  velocity  vector.  Note  that  Eq.  (21)  follows  simply  from 
conservation  of  momentum  and  energy  and  that,  to  this  point,  no 
approximation  has  been  made.  If  molecules  1  and  2  initially  hit  head-on, 
then  the  deflection  of  the  product  molecule  3  from  the  initial  center  of 
mass  velocity  direction  can  be  represented  by  the  angle  a,  where 


cos(a) 


v3‘vcm 

V3vcm 


If  cos(a)  is  computed  from  Eq.  (21),  there  results 


cos (a) 


vcm  +  CVrf^-1) _ 

V(v2m  ♦  e2v2f  ♦  2€VcmVrf(2C-l)) 


where  e  is  the  mass  ratio 


m3  +  m4 


and,  for  a  hard  sphere  collision,  t  is  a  random  variable  in  the  range  of  0 
to  1.  For  a  non  hard  sphere  collision  (the  real  case,  of  course),  Eq.  (23) 
still  applies,  but  t  is  no  longer  evenly  distributed  in  its  range.  In 
either  case,  Eq.  (23)  can  be  differentiated  with  respect  to  C  to  find  the 
maximum  possible  deflection  angle  (or  minimum  possible  cos(a) ) .  The  result 
is  the  relatively  simple  expression: 


(cos(a) ]min 


/I  -  €2(Vrf/Vcm)' 


Equation  (25)  is  in  a  form  which  can  be  easily  evaluated  for  the  casa  at 
hand,  if  molecules  1,  2,  3,  and  4  are  associated  with  0,  H2 ,  OH  and  H, 
respectively.  Note  that  e  is  small  (1/18),  which  means  that  the  maximum 
possible  deflection  angle  is  not  very  large.  This  is  consistent  with  the 
intuitive  reasoning  that  the  OH  molecule  cannot  deviate  far  from  the 
trajectory  of  the  original  0  molecule,  because  the  H  molecule  is  so  light 
and  cannot  compensate  greatly  in  the  center  of  mass  frame.  When 
appropriate  numbers  are  inserted  for  the  case  of  interest,  there  results: 
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where  nQ  is  the  atomic  oxygen  number  density.  Each  reaction  will  scatter 
back  and,  it  is  assumed,  spread  the  OH  products  evenly  over  a  cone  of 
half-angle  a^,,,.  The  cone  base  will  have  an  area  of  ir(L+s)  2tan2(amax) , 
where  L  is  a  characteristic  dimension  of  the  separation  of  the  orifice  from 
the  aperature.  It  is  assumed  that  L  is  on  the  order  of  1  meter.  The  total 
number  of  molecules  returned  per  unit  area  per  unit  time  from  this  region, 
AnQH,  can  be  expressed 


An, 


BdflAsn0<r  Vr 
0H  ir(L+s)2tan2(amax) 


(32) 


If  this  relation  is  integrated  from  s  equal  0  to  ®,  the  total  return  flux 
of  OH,  qQH,  is  evaluated  to  be 


qQH  2«  6.7xl07  molecules/cm2/ s 


(33) 


Hence,  if  a  monolayer  is  assumed  to  be  approximately  1015  molecules  per 
square  centimeter,  a  monolayer  will  be  formed  every  170  days  or, 
approximately  two  monolayers  a  year.  It  is  felt  that  this  level,  if 
accurate,  is  probably  small  enough  to  be  unimportant. 

2 . 5  Significant  Factors  Affecting  the  Return  Flux  Estimate 


Before  the  computed  flux  is  accepted,  and  the  effect  is  deemed 
unimportant,  it  is  helpful  to  consider  the  errors  in  Eq.  (33).  As  was 
mentioned  in  the  text,  the  reactive  cross  section  was  based  on  an 
extrapolation,  and  the  initial  H2  temperature  was  a  guess.  The  major 
source  of  error,  however,  is  the  atomic  oxygen  number  density,  nQ.  The 
value  used  was  4xl04  molecules/cm2,  corresponding  to  an  exospheric 
temperature  of  about  1000  K.  If  the  actual  exospheric  temperature  were  to 
be  2000  K  (corresponding  to  high  solar  activity),  the  atomic  oxygen  number 
density,  and  the  estimated  return  flux,  would  rise  by  more  than  two  orders 
of  magnitude.  Hence,  although  the  return  flux  is  computed  to  be 
negligible,  that  conclusion  could  change  if  the  exospheric  temperature  rose 
substantially  above  the  assumed  value  of  1000  K. 
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3.  APERATURE  AND  MIRROR  FLUXES  OF  ATMOSPHERIC  0  AND  H 


At  the  900  km  altitude  of  interest,  atomic  hydrogen  has  a  number 
density  comparable  to  that  of  atomic  oxygen.  (For  an  exospheric 
temperature  of  1000  K,  n^  2<  4xl04  molecules /cm"*  and  nH  ^  8xl04.)  Since  the 
atomic  oxygen  will  stick  to  cryogenic  surfaces,  and  the  fluxes  of  the  two 
species  are  comparable,  it  is  possible  that  surface  reactions  will  form  OH 
on  the  mirror.  Note  that  this  source  of  contamination  is  purely 
atmospheric,  and  does  not  involve  the  H2  coolant. 

A  crucial  element  of  the  present  calculation  is  the  determination  of 
what  fraction  of  the  flux  incident  on  the  aperature  actually  makes  it  to 
the  mirror.  For  a  variety  of  conditions,  two  fluxes  were  calculated:  a) 
the  flux  at  the  aperature,  and  b)  the  direct  flux  to  the  mirror.  These  two 
fluxes  bracket  the  possible  flux  to  -  the  mirror.  If  the  baffles  are  such 
that  all  molecules  impinging  on  them  stick,  then  the  direct  mirror  flux 
represents  the  total.  However,  if  a  portion  of  the  molecules  that 
initially  hit  the  baffles  bounce  around,  then  they  can  eventually  reach  the 
mirror. 

The  procedure  that  was  used  was  to  introduce  molecules  at  the 
aperature  appropriate  to  a  Maxwellian  velocity  distribution,  as  is 
described  in  Ref.  8.  The  trajectories  were  then  traced  to  see  what 
fraction  actually  impinged  on  the  mirror.  Cases  were  run  for  atomic  oxygen 
and  atomic  hydrogen;  for  exospheric  temperatures  of  500,  1000,  and  1500  K; 
and  for  varying  angle  between  the  telescope  centerline  and  the  oncoming  ram 
direction.  The  case  for  T^  ■  500  K  had  negligible  atomic  oxygen,  so  that 
is  not  shown  here.  Figures  1  and  2  show  the  O  and  H  flux,  respectively, 
for  T^  ■  1000  K,  and  Figures  3  and  4  show  the  same  quantities  for  T^  *  1500 
K. 

Figure  1  shows  that  when  the  telescope  is  facing  into  the  ram,  the 
atomic  oxygen  flux  at  the  mirror  is  reduced  by  less  than  50%  at  Ta  *  1000  K 
compared  to  the  flux  at  the  aperature.  Most  of  the  oxygen  atoms  go 
straight  through  to  hit  the  mirror.  As  the  telescope  points  away  from  ram, 
however,  the  flux  that  hits  the  mirror  drops  much  more  precipitously  than 
the  flux  at  the  aperature.  This  is  easy  to  understand  physically,  since 
the  mean  flow  velocity  becomes  pointed  away  from  the  direction  necessary  to 
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impact  the  mirror  directly.  The  mirror  flux  is  atill  about  109 
raolecules/cm2/s  at  20  degrees  from  ram,  which  is  quite  likely  still 
significant.  By  30  degrees  from  ram  the  mirror  flux  has  dropped  two  more 
orders  of  magnitude,  and  is  probably  no  longer  important. 

Figure  2  shows  the  H  flux  for  the  same  atmosphere.  In  the  ram 
direction  the  mirror  flux  is  an  order  of  magnitude  down  from  the  aperature 
flux,  which  is  a  much  larger  factor  than  was  observed  for  the  atomic  oxygen 
flux  above.  The  reason  for  this  is  the  much  larger  spread  in  the  velocity 
distribution  for  the  atomic  hydrogen.  This  same  spread,  however,  makes  the 
mirror  flux  drop  off  much  more  slowly  with  angle  from  ram  than  did  the 
atomic  oxygen.  There  was  even  calculable  flux  to  the  mirror  at  130  degrees 
from  ram,  though  it  is  not  significant. 

As  the  exospheric  temperature  is  raised  in  Fig3.  3  and  4,  changes  are 
observed  in  both  the  fluxes.  Firstly,  there  is  much  more  atomic  oxygen  at 
the  higher  temperature,  and  approximately  the  same  amount  of  atomic 
hydrogen.  Hence,  the  atomic  oxygen  curves  are  shifted  upwards  relative  to 


0  FLUX  FOR  T  =1000K 


ANGLE  FROM  RAM  (DEG) 


-  RT  APERTURE 

-  RT  MIRROR 


Figure  1.  The  Aperature  and  Direct  Mirror  Flux  of  Atomic  Oxygen 
Calculated  for  an  Exospheric  Temperature  of  1000  K  as  a 
Function  of  Angle  from  Ram. 
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H  FLUX  FOR  T  -1000K 
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Figure  2.  The  Aperature  and  Direct  Mirror  Flux  of  Atomic  Hydrogen 
Calculated  for  an  Exospheric  Temperature  of  1000  K  as  a 
Function  of  Angle  From  Ram. 


0  FLUX  FOR  T  -1500K 
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Figure  3.  The  Aperature  and  Direct  Mirror  Flux  of  Atomic  Oxygen 
Calculated  for  an  Exospheric  Temperature  of  1500  K  as  a 
Function  of  Angle  From  Ram. 


11 


H  FLUX  FOR  T  =1500K 


ANGLE  FROM  RAM  (DEG) 


FIT  RPERHTURE 
RT  MIRROR 


Figure  4.  The  Aperature  and  Direct  Mirror  Flux  of  Atomic  Hydrogen 
Calculated  for  an  Exospheric  Temperature  of  1500  X  as  a 
Function  of  Angle  From  Ram. 


the  atomic  hydrogen.  Also,  the  increased  spread  in  the  velocity 
distribution  function  caused  by  the  higher  temperature  produces  a  somewhat 
larger  drop  between  aperature  and  mirror  fluxes  at  ram  for  both  species. 
For  Tgj  *  1500  X,  the  telescope  has  to  be  approximately  forty  degrees  from 
ram  before  the  O  flux  drops  to  107,  compared  to  thirty  degrees  at  the  lower 
temperature . 

These  results  again  point  out  the  great  sensitivity  to  atmospheric 
temperature,  but  it  can  also  generally  be  concluded  that  facing  into  the 
ram  for  long  periods  of  time  is  likely  to  be  a  hazardous  operation.  If  it 
is  done  at  all,  it  should  probably  be  at  the  end  of  the  mission. 
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